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BRIEF  ABSTRACT 


One  of  the  prime  difficulties  in  developing  two  dimensional  dynamic 
elastic  plate  theories  from  the  three  dimensional  equations  of  elasticity  is 
the  choice  of  functional  dependence  on  the  thickness  coordinate. 

This  difficulty  may  be  circumvented  by  formulating  the  problem  first  as  a 
boundary  Integral  equation;  then  the  dependence  on  the  independent  variable 
through  the  plate  thickness  follows  form  as  a  direct  quadrature  with  no 
assumptions  of  functional  form  required. 

In  particular,  the  examination  of  separate  symmetric  and  antisymmetric 
modes  allows  the  boundary  integral  equation  to  be  written  with  unknowns 


evaluated  on  a  single  surface. 
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Consider  an  infinite  elastic  plate  of  thickness  h.  The  dispersion 
equation  for  elastic  vibrations  of  such  a  plate  is  well  known,  e.g.  Ewing, 
Jardetsky  and  Press  (1957).  The  first  step  in  this  research  is  to  duplicate 
such  a  dispersion  equation  by  means  of  a  boundary  integral  equation  approach. 
The  purpose  of  this  development  is  not  to  demonstrate  an  alternative 
derivation  procedure  (in  fact,  this  approach  would  appear  clumsy  compared  to 
the  classical  one) ,  but  rather  to  illustrate  an  approach  which  may  have  much 
wider  application.  In  particular,  the  boundary  Integral  equation  method 
allows  solution  of  the  exact,  three  dimensional  equations  of  linear  elasticity, 
but  on  two  dimensional  surfaces.  This  accomplishes  the  primary  goal  of  plate 
theory,  i.e.  the  reduction  of  the  dimensionality  of  the  governing  equations, 
without  paying  the  corresponding  price  of  a  physical  approximation. 

Furthermore,  if  approximations  are  made  based  on  the  physical  property  that 
the  plate  thickness  is  small  compared  to  the  wave  lengths  involved,  their 
effect  can  be  followed  directly  in  the  governing  integral  equations  which 
may  themselves  be  solved  through  an  expansion  in  orders  of  plate  thickness. 

Define  the  upper  and  lower  surfaces,  at  y  •  +h/2  and  -h/2 
respectively,  as  Sy  and  S^.  These  shall  be  taken  to  be  stress  free. 

The  problem  may  be  treated  as  two  dimensional  in  that  waves  are  taken  to 
propagate  only  in  one  horizontal  direction  (x)  with  no  dependence  on  the 
second  horizontal  direction  (z  normal  to  x  and  y) .  This  allows  the 
introduction  of  two  scalar  displacement  potentials,  p  and  p  ,  e.g. 


u  ■  Vp  +  Vx(pk) 

u  ■  p,  +  p, 
x  T*x  Yfy 
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These  displacement  potentials  both  satisfy  wave  equations.  Assuming 
time  harmonic  dependence  exp  (-iwt) ,  these  equations  are  then  Helmholtz 
equations » 

V2*  +  k2  *  -  0 
V2*  +  k2  *  -  0 

where  the  dilatational  and  transverse  wave  numbers  ,  and  k^,  are 
w//(A  +  2p)/p  and  w//y/p  respectively  with  Lamd  constants  A  and  p. 

The  stresses  are  related  to  the  displacement  potentials  by 

°xx  *  X[*'xx  +  ♦•yyl  +  ^’xx  +  *»xyl 
°yy  *  ^^»xx  +  ^*yyl  +  2yl9»yy  ~  '('•xyJ 

°xy  *  w^*xy  +  ^»yy  "  ^»xx^ 


Boundary  conditions  for  stress  free  surfaces  at  constant  y  require  o 
and  a  to  be  zero  or 

yy 


xy 


*’*y  '  *•*«  +  kT  */2 
*>„  ■  ■*?  *'2  -  ♦>« 


The  governing  two  dimensional  Helmholtz  equation  may  be  rewritten  as 
an  Integral  (Weber)  equation,  e.g.  Sneddon  (1957)  by  use  of  Green's  theorem 
and  a  fundamental  solution.  (Alternative  formulations  are  available  which 
treat  stress  directly  rather  than  through  a  displacement  potential,  e.g. 
Cruse  and  Rizzo  (1968)).  For  an  infinite  elastic  plate,  the  Weber  equations 


are 


-  A  - 


where  S  represents  both  the  upper  and  lower  surfaces.,  Sy(y  ■  +h/2) , 

S  (y  =  -h/2).  n  is  the  outward  (from  the  solid)  normal  to  these  surfaces, 
li  o 

r  is  a  field  point,  rQ  is  an  integration  variable,  R  is  the  distance  between 
r  and  rQ  and  e  is  (0,  1/2,  1)  depending  on  whether  r  is  exterior  to,  on  the 
surface  of,  or  interior  to  the  elastic  body  respectively.  When  r  is  allowed 
to  approach  either  surface,  these  equations  reduce  the  dimensionality  of  the 
problem,  i.e.  they  Involve  only  surface  values  of  the  dependent  variables. 

The  variation  of  these  variables  through  the  thickness  of  the  plate  is 
prescribed  by  using  the  same  equations  with  the  field  point  interior  to  the 
plate;  here,  however,  once  the  surface  values  have  been  obtained  by  their  own 
system  of  equations,  the  interior  values  are  found  through  a  direct 
quadrature. 

Two  basic  forms  of  motion  are  possible;  the  symmetric  form  in  which 
~  ancl  t*ie  antisymmetric  form  in  which  ♦y  5  and  ♦y  = 
where  U,  L  refer  to  upper  and  lower  surfaces  respectively.  A  periodic 
dependence  on  x  in  the  form  exp.  (i*cx)  will  be  assumed;  this  form  could  also 
be  obtained  through  the  integral  equations. 

Consider  the  field  point  to  lie  on  the  upper  surface.  By  symmetry, 
only  $  and  ^  values  on  this  surface  are  required  for  the  complete  problem. 
Define 

$(x,h/2)  -  ♦y(x)  "  ±4>L(x)  ■  A  exp  (iicx) 
f(x,h/2)  ■  il'y(x)  -  T^l(x)  -  B  exp  (iicx) 

3*  3Vx)  *L 

(x,h/2)  -  — — —  -  +  (x)  -  C  exp  (iicx) 

d<b 

(x,h/2)  -  ~  (x)  -  ±  7~  (x)  -  D  exp  (iicx) 
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where  the  two  signs  refer  to  symmetric  and  antisymmetric  modes  respectively, 
and  the  two  subscripts  U  and  L  refer  to  uppper  and  lower  surfaces 
respectively. 

The  equation  then  becomes 

3*, 


<l/2)+„(»)  -  (1/41)  £_  (*o)]d  *o 

+  <l/4i)  L 

*  ♦o(,o)  a^["o1,(kDRi)]] 


yo-  y  -  h/2 


d  *o 

'y  ■  h/2  ■  -y 


with  a  similar  equation  on  4*  with  k,j,  replacing  kp  and  the  signs  on  the  lower 
surface  Integral  reversed.  Here  Rp  *  Jx  -  xq|  ■  Jc(»  ■  |Jx  -  x0j2  +  h2J^2 
and  C  •  x0  -  x;  where  9Rp^/3yo  vanishes  on  the  upper  surface.. 

The  boundary  conditions  provide 

C  -  iicB  -  |i  k*/2ic}B 
D  -  -ixA  +  ^i  k2^2tcjA 


Then 


end 


(1/2) A  exp 


+(l/4i>  f 


(Ikx)  -  (1/41) J  exp  [ixxjH^^kpRpjjd  x0 

y  *  y0 

•A(V/l4.)Hi1)(kDlV)  +  c  "ol>KRi.)|*,p(1"1.)d  *1 

'y  -  h/2  -  -y 


h/2 


(1/2)B  exp (lxx)  -  (l/4i)j  l-D  exp  JixxjH^  Id  x 


'  I 


y  •  y  -  h/2 

O 


ta/4I,LFB-(kT5A.)"il>(kTi<i.)  +  D  "o1>(kA)l“'(1".)d 


J 

Jy 


-  h/2  •  -y 


These  can  be  shifted  to  F,  =  xq  -  x,  thus  removing  the  x  dependence  from  the 
left  hand  side  of  the  equation.  These  are  clearly  Fourier  Transforms  and  may 
be  evaluated  using  the  Bateman  tables,  Erdelyi,  et  al.  (1954),  e.g*: 


,00 

expdicOH^^kUljdS 

*  —00 

-  2[  coadCJH^OcOdS 
'  0 

(k2  -  K2)-l/2.  0  <  k  <  k 

-i(<2  -  k2)-1^2;  k  <  <  <  ®j 


)  exp(iieC)Ht:i)(k[C2  +  h2]l/2)d£ 

J-co  ° 

fCO 

-  2 [  cos(icC)Ha)(k[C2  +  h2]l/2)de 

J  o  ° 

|+  (k2  -  <2)“^2  exp  [ih(k2  -  k2)i/2];  0  <  k  <  k 
[-i(ic2  -  k2)"1/2  exp[-h(<2  -  k2)^2J;  k  <  ic  <  «• 

(-kh)exp(iK5)Hj1)(k[C2  +  h2]l/2)*(€2  +  h2)'1/2  d£ 

»  2^ (-khJcosdOH^dU2  +  h2]1/2)*^2  +  h2)"l/2  dC 
*  0 

+i  exp  [ih(k2  -  K2)1/2];  0  <  k  <  k 

-  2  ' 

+i  exp  [-h(ie2  -  k2)1'2];  k  <  ic  <  • 
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There  are  four  homogeneous  equations  on  the  four  unknown 
coefficients  which  then  must  define  a  dispersion  relationship. 


(1/2)A  -  (l/4i)  -2C 


+  02' 


N  - - ‘2H 

-if*2  -  k2)-1'2]  l+exp^-hp  -  k2)!/2] 


(‘2  -  »*) 


^k2  -  tc2j_1/2  exp|ih^kp  -  K2j1^  0 

-if*2  -  k2]-1/2  exp[-hp  -  k^j  l/2J  „ 


0  <  K  < 


<  K  <  » 


M‘2  -  “TT 

^k2  -  ie2j-1/2  exp£lh|k2  -  ic2^  l/2J  0  <  k  <  k^, 

-i|<2  -  k£j-1^2  exp£-h|tc2  -  k^j  1^2J  kj  <  k  <  • 


(1/2) B  -  (1/41)  -2D 


(k2  -  K2)-^2 

-i(K2-k2)- 


+B»2i- 


C-  [iK  -  ik2/2^B 
D--|iK  -  ikJ^icjA 

Since  <  k^,  there  are  three  distinct  regions  of  solution: 

I  :  0  <  k  <  kp  < 

II  :  0  <  kjj  <  k  <  kT 

III:  0  <  kjj  <  <  k 

in  addition  to  the  two  basic  modes  (symmetric  and  antisymmetric) .  Only  case 
(ill)  will  be  considered  since  only  long  waves  are  anticipated.  Consider  the 
case  (III)  -  symmetric:  the  governing  equations  lead  directly  to  the  well- 
known  dispersion  equation,  for  general  elastic,  e.g.  Ewing,  Jardetsky  and 
Press  6-12  (1957),  with  v2  -  ic2  -  k2  and  v’2  -  k2  -  k2, 


Avv  *  K2 


tanh  [hvt/2]  _ 

tanh  [hv/2]  "  (ic2  +  v'2)2 


For  h  very  small ,  this  reduces  to 

A(!  -  k2/k2)(K2/k2)  -  1 


which  is  the  well-known  elastic  plate  dispersion  equation,  e.g.  Ewing, 

Jardetsky  and  Press,  6-16,  (1957). 

Case  III  -  antisymmetric  follows  in  much  the  same  way,  resulting  in  the 

dispersion  equation 

tanh  [hv/2]  _  4vvf  x2 

tanh  [hv'/2]  (ic2  +  v’2)2 


as  found  as  equation  6-11  In  Ewing,  Jardetsky  and  Press  (1957)* 

Interior  values  of  $  and  ip  can  be  found  once  the  surface  values  are 
known  (at  least  to  within  a  constant  multiplier).  Placing  the  field  point  at 
(x,y)  where  -h/2  <  y  <  +h/2  yields 


♦<*,y>  -  (l/4i)  ♦„(*„) "‘“(Wo) 


y„  - 


d  x 

o 


y 


o 


«  -h/2 


and  a  similar  equation  on  t|».  Using  slight  modifications  of  the  previous 
Fourier  transforms  leads  directly  to  interior  values,  e.g.  the  symmetric 
solution  for  4  for  case  III  is 
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♦(x,y)  ■  -exp(itcx)  •  exp^h^K2  -  lc2^2  2J 


•  [*vp  ■  •  Wr*  -  ty'i 


which  again  agrees  with  equation  6-13  of  Ewing,  Jardetsky  and  Press. 
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